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With the advances in miniaturization, understanding and controlling properties of significant 
technological systems like silicon in nano regime assumes considerable importance. The small silicon 
clusters in the size range of 15-20 atoms are known to fragment upon heating. In the present work 
we demonstrate that it is possible to stabilize such clusters by introducing appropriate dopant (in 
this case Ti). Specifically, by using the first principle density functional simulations we show that 
Ti doped Sii6, having the Frank-Kasper geometry, remains stable till 2200 K and fragments only 
above 2600 K. The observed melting transition is a two step process. The first step is initiated 
by the surface melting around 600 K. The second step is the destruction of the cage which occurs 
around 2250 K giving rise to a peak in the heat capacity curve. 

PACS numbers: 31.15.Ew, 31.15.Qg, 36.40.Ei, 36.40.Qv 



I. INTRODUCTION 

With the advent of nano science and technology, in- 
vestigating physical properties of clusters has assumed 
considerable importance^ As the size of the system be- 
comes smaller and smaller its physical and chemical prop- 
erties are expected to change, sometimes significantly. In 
fact, this is what drives the quest for nano materials re- 
search. However, it may happen that this size effect could 
drive a property towards an undesirable direction. The fi- 
nite temperature properties of small silicon clusters is one 
such example. In our recent finite temperature studies^ 
it was found that the small clusters of silicon (Si„, n=15 
and 20) become unstable and fragment when heated up to 
about 1500 K. Typically, our extensive density functional 
simulations showed that Siis fragments around 1800 K 
in to Sig and Sig, although we also observe Siio+Si5 and 
Si8+Si7 for a short time span. It exists in a liquid like 
phase over a short temperature range (900 to 1400 K) 
before fragmenting. On the other hand, Si2o transforms 
from the ground state to the first isomer (with two dis- 
tinct Siio units) around 1000 K and eventually fragments 
in to two Siio units around 1200 K. Thus, Si2o does not 
show any solid-like to liquid-like transition prior to frag- 
mentation. The result although not surprising, can have 
technological implications, silicon being the key ingredi- 
ent in most of the semiconductor devices. The present 
work addresses the issue of remedying this difficulty by 
using appropriate dopant. That the impurities induce 
significant changes in the structure as well as electronic 
properties is well known in the solid state community. 
In the context of clusters several studies have been re- 
ported. For example, single atom of tin or aluminum in 
lithium clusters are known to change the nature of bond- 
ing among the host atoms In a recent work Mottet et. 
al& demonstrated that a single impurity (Ni or Cu) in 
Ag55 cluster changes structural and thermal properties 
dramatically. Interestingly, this effect was found to be 
significant only on an icosahedral host. 



Silicon is the most important semiconducting mate- 
rials in the microelectronics industry. In medium-sized 
regime, the structure and properties of materials often 
differ dramatically from those of the bulk. Over the past 
several years there have been extensive theoretical and 
experimental work describing the ground state geome- 
tries of medium sized silicon clusters Ai&SiiSiii These 
clusters are built upon stable tricapped trigonal prism 
(TTP) units which plays a crucial role in the finite tem- 
perature behavior of the clusteriii The high stability of 
TTP units results in to fragmentation of these clusters 
instead of solid-like to liquid-like transition £ The ground 
state geometries of these clusters are well established in 
the literature. The ground state geometry of Sii6 can be 
described as two fused pentagonal prisms. 8 

Experiments on doped silicon clusters Si„M ( M=Ti, 
Hf, Cr, Mo, and W) have shown large abundances for n — 
15, 16 and low intensities for other sizesii2ii2ii4 Specifi- 
cally the abundances drop drastically beyond n = 16, 
giving support to the predictions of the exceptional sta- 
bility of n = 16 clusters. Vijay Kumar and co-workersi& 
showed that for n = 8 — 12 basket-like open structures 
to be the most favorable, while for n — 13 — 16, the 
metal atom is completely surrounded by silicon atoms. 
Recently, measurements of electron affinities of Ti-doped 
silicon clusters** have been found to have a minimum 
at n = 16, suggesting the closed electronic shell nature 
and strong stability of Si^Ti. The electronic structure 
of silicon clusters in the size range noted above has been 
studied by ab initio methodsAiS It may be noted that 
none of the silicon clusters in this size range form stable 
caged structures. In a quest to obtain stable silicon cage 
structures, Vijay Kumar and co-workersi^iSiiii^ found 
that it is possible to stabilize a caged structure for clus- 
ter using a class of dopants. Specifically, they showed 
that a single impurity of transition metal atoms like Ti, 
Zr, Hf enhances the binding energy of Sii6 and changes 
the structure to a caged one, similar to carbon cages. 
The gain in binding energy is substantial and is of the 
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order of few eV and HOMO-LUMO gap is more than 
1 eV which make these structures stable. Among the 
systems studied, Sii 6 Ti has the largest HOMO-LUMO 
gap (2.35 eV) as compared to other dopants of Sii6 like 
Zr and Hfi^ The ground state geometry for Si^Ti has 
been found to be a Frank Kasper polyhedron, which is 
also shown in Fig. Q}a. Therefore Si^Ti is a likely can- 
didate to show a stable behavior at finite temperatures. 
As we shall demonstrate, by adding a dopant like Ti, we 
could avoid the fragmentation observed in small Si clus- 
ters. The cluster shows a surface melting of Si atoms. 
However, the motion of Si atoms is restricted to a small 
shell around Ti until 2200 K. The cluster is on the verge 
of fragmentation around 2600 K, almost 1000 K higher 
than that in pure Si clusters. 



II. COMPUTATIONAL DETAILS 

We have carried out isokinetic Born-Oppenheimer 
Molecular Dynamic (BOMD) simulations using ultra- 
soft pseudopotentials within the Generalized Gradient 
Approximation (GGA), as implemented in the VASP 
packaged For computing heat capacities, the BOMD 
calculations were carried out for 17 different tempera- 
tures in the range of 100 K to 3000 K, each with the 
duration of 150 ps or more, which result in to a total 
simulation time of 2.6 ns. In order to get converged heat 
capacity curve especially in the region of coexistence, 
more temperatures were required with longer simulation 
times. We have discarded the first 30 ps of each tem- 
perature for thermalization. To analyze the thermody- 
namic properties, we first calculate the ionic specific heat 
by using the Multiple Histogram (MH) technique^fiiSl 
We extract the classical ionic density of states (ft(E)) 
of the system, or equivalently the classical ionic entropy, 
S(E) = k B \nQ(E) following the MH technique. With 
S(E) in hand, one can evaluate thermodynamic averages 
in a variety of ensembles. We focus in this work on the 
ionic specific heat and the caloric curve. In the canon- 
ical ensemble, the specific heat is defined as usual by 
C(T) = 8U(T)/dT, where U(T) = J Ep(E,T) dE is the 
average total energy, and where the probability of ob- 
serving an energy E at a temperature T is given by the 
Gibbs distribution p(E, T) = fl(E) exp{-E/k B T)/Z{T), 
with Z(T) the normalizing canonical partition function. 
We normalize the calculated canonical specific heat by 
the zero-temperature classical limit of the rotational plus 
vibrational specific heat, i.e., Cq — (3N — 9/2)fce. 

We have calculated the root-mean-square bond length 
fluctuations (<5 rms ) for Si-Si and Si-Ti bonds separately to 
compare the difference between them. The <5 rms is defined 
as 
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where N is the number of bonds in the system(N=120 
for Si-Si, and N=16 for Si-Ti), is the distance between 



atoms i and j, and (...)* denotes a time average over the 
entire trajectory. Mean Square Displacement (MSD) is 
another traditional parameter used for determining phase 
transition and is defined as, 
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here N is the number of atoms in the system(N=16 for 
Si and N=l for Ti) and R is the position of the I th atom. 
Here we average over M different time origins to m span- 
ning over the entire trajectory. The interval between the 
consecutive to m for the average was taken to be about 
0.3 ps. The MSD of a cluster indicate the displacement 
of atoms in the cluster as a function of time. 

We have also calculated radial distribution function 
(g(r)) and deformation parameter (sdef)- g( r ) is defined 
as the average number of atoms within the region r and 
r+dr. The shape deformation parameter (£def) is defined 
as, 
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where Q x > Q y > Q z are the eigenvalues, in descending 
order, of the quadrupole tensor 



Qij — Rn Rij 
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Here i and j run from 1 to 3, I runs over the number of 
ions, and Rn is the \ th coordinate of ion / relative to the 
COM of the cluster. A spherical system((5 a; — Q y = Q z ) 
has ed e /=l and larger values of e e i e f indicates deviation 
of the shape of the cluster from sphericity. We use the 
electron localization function (ELF) — to investigate the 
nature of bonding. For a single determinantal wave func- 
tion built from Kohn-Sham orbitals ipt, the ELF is de- 
fined as 
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with p = p(r) is the valence-electron density. The ELF is 
defined in such a way that its value is unity for completely 
localized systems and 0.5 for homogeneous electron gas. 



III. RESULTS AND DISCUSSION 

We begin our discussion by noting the ground state 
(GS) and some low lying structures of Si^Ti which are 
shown in Fig. ^ As has been already mentioned the 
ground state of Si^Ti is Frank-Kasper polyhedron. The 
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FIG. 1: The ground state geometry and few interesting iso- 
mers of Sii6Ti 



central Ti atom is surrounded by 16 Si atoms within two 
closely spaced shells one with 12 atoms, all equidistant 
from central Ti atom and another shell consisting of 4 Si 
atoms forming a tetrahedron with Ti at the center. The 
first isomer (Fig. GJb) is energetically very close to the 
ground state where all Si atoms are nearly equidistant 
from central Ti atom (A E = .009 eV). Isomers with ei- 
ther distorted cage of Si atoms or without Si cage are 
quiet high in energy compared to the lowest two struc- 
tures. Isomer shown in Fig. ^f occurs at very high tem- 
perature (around 2600 K) and indicates possible path for 
fragmentation. At this point its interesting to compare 
pure silicon clusters with doped clusters. Si clusters have 
very interesting growth pattern. The ground state ge- 
ometries of Si n till n = 22 are prolate and for larger sizes 
the GS transforms into spherical structures. It has been 
also observed that a TTP unit is a part of prolate ground 
states. Presence of highly stable TTP unit also bypass 





FIG. 2: The isosurfaces of ELF at Xelf=-70 and Xelf=0.55 
are shown in (a) and (b) respectively. 
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FIG. 3: The heat capacity for Si^Ti computed over last 120 
ps. The peak is at 2250 K. 



the liquid-like state and leads to fragmentation of pure 
Si„ clusters around 1200 K to 1800 K in the size range 
of 15-20. Further, the ground state of Sii6 is a prolate 
structure and a cage like isomer is found to be difficult 
to stabilize. 

The nature of bonding in SiieTi can be discussed by 
examining the iso-surfaces of total charge density and 
ELF. In Fig. |21 we show the iso-surface of ELF for value 
of 0.7 (Fig. Ha) and 0.55 (Fig.^b). At higher values of 
Xelf the figure clearly shows localization of the charge on 
hexagonal rings depicting the covalent nature of bonding. 
These atoms are connected to each other by the shortest 
bonds (2.37-2.43 A) and are connected to the remain- 
ing four Si atoms at very low value of xelf ~ 0.60 ( 
bondlcngths around 2.65 A). The bonding between sili- 
con atoms belonging to two different shells is metallic like 
in the sense it is due to delocalized charge distribution. 

Next we show the ionic heat capacity of SiieTi in Fig. 
The heat capacity displays a broad melting peak at 
2250 K characteristic of a finite size system. A more 
careful observation of ionic motions and traditional pa- 
rameters like <5 rms reveals that the main peak in the heat 
capacity is related to complete breakdown of Si cage and 



500 1000 1500 2000 2500 3000 

Temperature (K) 



30 



60 

Time (ps) 



90 



120 



FIG. 4: The 5 Iula for Si-Si bonds (online red color) and Si- 
Ti bonds (online blue color) are shown. The region between 
600 K to 1800 K (shown by arrows for Si-Si bonds) corre- 
sponds to 'restricted' liquid-like behavior of Si cage. The 
rise in <5 rma observed after 2000 K (shown by arrow for Si-Ti 
bonds) is associated with the breaking of Si cage and diffusion 
of Ti through out the cluster. 



escape of Ti atom from the cage. The S lms for Si-Si and 
Si-Ti bonds, shown in Fig. 0] throws light on the finite 
temperature behavior of this cluster. The value of S lms 
for Si-Si bonds is considerably high at much lower tem- 
peratures (around 200 K) whereas the value of <5 rms for 
Si-Ti bonds is quite low till 1800 K. The higher value of 
<5rms for Si atoms indicates that Si atoms are mobile at 
low temperatures. More detailed analysis of ionic motion 
reveals that the cluster undergoes isomerization around 
as low temperature as 100 K. The cluster transforms from 
the GS to the first isomer and back to the GS which in 
turn results in to diffusion of at least 4 atoms. The MSD 
also support this observation. In Fig. and Fig. the 
MSD computed for Si and Ti atoms at some relevant tem- 
peratures are shown respectively. As can be noted from 
Fig. El the MSD for Si atoms saturate around 600 K indi- 
cating that the atoms are diffusing on the spherical shell 
around Ti atom and the cluster is partially melted. On 
the contrary, as can be seen from Fig.|S|the Ti atom shows 
very small displacement till 2000 K. Although the MSD 
shows liquid-like behavior around much lower tempera- 
ture (600 K) the Si atoms are confined to a small shell 
around Ti (preserving the spherical shape) until quite 
high temperatures. Its only around 2000 K or so that 
the Si cage breaks down giving rise to a peak in the heat 
capacity around 2200 K. That the motion of Si atoms 
is restricted to a shell is clear from the examination of 
radial distribution function. 

In Fig. [7] we show the radial distribution function cal- 
culated from the center of mass (Ti site) for some repre- 
sentative temperatures. The peak near origin is due to 
the Ti atom whereas a single peak around 2.8 A is due 
to Si cage. The width of this peak does not change sig- 
nificantly till 1600 K. A continuous distribution emerges 
after 2200 K indicating that the Si cage is destroyed and 



FIG. 5: The MSD calculated for Si atoms for few relevant 
temperatures. 
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FIG. 6: The MSD calculated for Ti atom for five different 
temperature 



Ti has escaped from the cage. This change is accompa- 
nied by the shape change of the cluster. 

In Fig. |21 we show shape deformation parameter as a 
function of temperature. It can be seen that till 2000 K 
value of £def is nearly 1 indicating a spherical shape with 
restricted motion of Si atoms on the surface. The Edef 
increases significantly above 2000 K indicating the break- 
down of Si cage. Our simulation shows that a fragmenta- 
tion channel opens up around 2600 K possibly because of 
emergence of central Ti atom on the surface. We have ob- 
served that Si^Ti around this temperature fragments in 
to two parts, the most probable fragments are Sii2Ti+Si4 
and Sii3Ti+Si3. A typical configuration leading to frag- 
mentation of SiigTi is also shown in the Fig. [SJ 



IV. SUMMARY AND CONCLUSION 

To summarize, our DFT simulations clearly demon- 
strate that its possible to avoid fragmentation of small 
Si clusters by doping them with appropriate impurity, in 
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FIG. 7: The radial distribution function calculated for five 
different temperatures. 



this case Ti. Contrary to the pure Si clusters, the doped 
cluster undergoes a solid-like to liquid-like transition, and 
remain stable at least up to 2600 K. However, the 'melt- 
ing' occur in two steps, the first one is a low temperature 
process (around 600 K), where Si atoms diffuse on the 
shell, whereas the peak in the heat capacity (2250 K) is 
associated with the breaking of Si cage and escape of Ti 
atom from the cage. The beginning of the fragmentation 
is observed around 2600 K which is almost 1000 K higher 
than that observed in pure Si clusters in this size range. 
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FIG. 8: The deformation parameter is plotted as a function 
of temperature. Note the sudden rise in the deformation pa- 
rameter after 2000 K. A typical configuration shown in the 
Fig occurs around 2400 K, and indicates possible path for 
fragmentation. 
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